Phylogenetic tree reconstruction from genome-scale 
metabolic models 

D. Gamermann* 1,2 , A. Montagud 2 , J. A. Conejero 2 , P. F. de Cordoba 2 and J. F. Urchieguia 2 

1 Catedra Energesis de Tecnologi'a Interdisciplinar, Universidad Catolica de Valencia San Vicente Martir, 
Guillem de Castro 94, E-46003, Valencia, Spain. 

2 lnstituto Universitario de Matematica Pura y Aplicada, Universidad Politecnica de Valencia, 
Camino de Vera 14, 46022 Valencia, Spain. 

Email: Daniel Gamermann*- daniel.gamermann@ucv.es; 

"Corresponding author 

Abstract 

A wide range of applications and research has been done with genome-scale metabolic models. In this work 
we describe a methodology for comparing metabolic networks constructed from genome-scale metabolic models 
and how to apply this comparison in order to infer evolutionary distances between different organisms. Our 
methodology allows a quantification of the metabolic differences between different species from a broad range of 
families and even kingdoms. This quantification is then applied in order to reconstruct phylogenetic trees for sets 
of various organisms. 

1 Introduction 

Metabolic models at the genome-scale are one of the pre-requisites for obtaining insight into the operation 
and regulation of metabolism as a whole [l}j4] . Uses of metabolic models embrace all aspects of biotechnology: 
from food [5] to pharmaceutical [6] and biofuels [7|j9]- Genome-scale metabolic network reconstruction is, in 
essence, a systematic assembly and organization of all reactions which build up the metabolism of a given 
organism. It usually starts with genome sequences in order to identify reactions and network topology. This 
methodology also offers an opportunity to systematically analyse omics datasets in the context of cellular 
metabolic phenotypes. 
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Reconstructions have now been built for a wide variety of organisms, and have been used toward five 



major ends 10 : contextualization of high-throughput data [4,8,9 11 , guidance of metabolic engineering 



12 , directing hypothesis-driven discovery 17 , interrogation of multi-species relationships 18 , and network 



property discovery 15 



Nowadays, phylogenies are becoming increasingly popular, being used in almost every branch of biol- 



ogy 16 . Beyond representing the relationships among species in the tree of life, phylogenies are used to 



describe relationships between paralogues in a gene family [17 , histories of populations 18 , the evolutionary 



and epidemiological dynamics of pathogens [19 20 , the genealogical relationship of somatic cells during dif- 



ferentiation and cancer development 21 and even the evolution of language 22 . More recently, molecular 



phylogenetics has become an indispensible tool for genome comparisons 23 - 25 



A phylogeny is a tree containing vertex that are connected by branches. Each branch represents the 
persistence of a genetic lineage through time, and each vertex represents the birth of a new lineage. If 
the tree represents the relationships among a group of species, then the vertex represent speciation events. 
Phylogenetic trees are not directly observed and are instead inferred from sequence or other data. Phy- 
logeny reconstruction methods arc either distance-based or character-based. In distance matrix methods, 
the distance between every pair of sequences is calculated, and the resulting distance matrix is used for tree 



reconstruction. For a very instructive review, please refer to 16 



This work is organized as follows. In the next section we briefly explain the genome-scale models we 
work with and how we define a parameter for comparing two models. Following that section, we explain how 
we recover the phylogenetic tree from the comparison matrix obtained for many metabolic models. This 
will be done taking into account the minimum spanning tree of a non-directed connected weighted network 
associated to these metabolic models. In the section after that, we present the results and a brief study of the 
sensibility of the comparison parameter, and in the last section we present a brief summary and overview. 



2 Comparison between metabolic models 

In a recent paper [26] , it has been presented a method for automatically generate genome-scale metabolic 



models from data contained in the KEGG 27 database. The method consists in searching the database for 



genes and pathways present in an organism and download the corresponding set of chemical reactions. The 
algorithm filters isosenzymes, or other repeated reactions, and may add missing reactions to a given pathway 
using a probabilistic criterion based on the comparison of the organism's pathway with the same pathway 
in other organisms. In this work we are using data obtained from this platform, but the method described 
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can, in principle, be used with any set of metabolic models given that the compounds names in the models 
follow the same standard (the same compound has the same name in all models). 

The first step in our work is to construct for every metabolic model A a non-directed connected network 
Na — (Va,Ea) from the information contained in it. Here, Va stands for the set of vertex (nodes) of 
A and Ea for its set of edges. A metabolic model comprises a set of chemical reactions. Each chemical 
reaction associates a set of substrates with a set of products. For constructing the network, first we define 
the set of vertex Va as the set of compounds in A (metabolites present in the model), assigning a vertex to 
each metabolite. The chemical reactions in the model will define the edges (links) of the network. If two 
metabolites appear as substrate and as a product, respectively, in a chemical reaction, an edge connecting 
the correspondent vertex is added to the network. A typical metabolic model of a prokaryote, with around 
1000 metabolites and the same number of chemical reactions, becomes through this process a non-directed 
connected network with 1000 vertex and around 3000 links. 

The problem at hand now is to elaborate a method for systematically compare and quantify the differences 
between two metabolic networks. For this purpose we define a parameter that scales between zero and infinity, 
zero meaning identical networks and infinity for networks that either share no node or no edge in common. 
The definition of this parameter is based on the identity of the nodes (the compounds), but not directly on 
the chemical reactions of the metabolic models, only indirectly through the edges of the network. 

Here we start with the metabolic networks of two organisms A — (Va,E a ) and B = (V b ,Eb)- The set 
of all metabolites in between the two organisms A U B = (Va U Vb, Ea U Eb) can be divided into a partition 
of three disjoint sets: the set of metabolites only present in A, the set of metabolites only present in B, and 
the set of metabolites common to both organisms: 



Vaub = (v A \v B )u(v A nv B )u(v B \v A ) (1) 



Only in A Common Only in B 



where \ stands for the difference of sets. A representation of this situation is shown in Figure 5.1 As it is 
represented there, each metabolite may have connections to metabolites within its set and connections to 
metabolites in the other sets. 

Suppose that Va U V b = {v±, . . . , v n }. Fix an arbitrary node Vi, 1 < i < n. We can consider its degree 
in A U B, i.e. the total number of connections of Vi to the rest of metabolites of Va U V b , that we denote 
by deg(vi). We can also consider the degree of vi when we restrict ourselves to the subnetwork generated 
by the vertex in (Va \ Vb) U {Uj}, that we will call deg A \ B (vi). Similarly, we can also define deg AnB (vi), 



and dcg B \ A (vi). With these degrees we can define, for each metabolite Uj € Va U Vb, the rate pa\ba of 
connections of to metabolites inside ^4 and not in B respect to the total number of connections of Vi, that 
is: 

degA\s(«i) 
deg^) 

Analogously, we can define 

degg^i) , deg AnB (u i ) 
deg(Vi) deg{vi) 

The following weighted sum of the rates Pa\b,i provides a parameter of the differentiation of A U B 
respect to A: 

\ V ^V B \ E A ^(v 3 )\ E degg) 

On the one hand, the rates Pa\b,i are multiplied by the inverse of the total number of connections of Vi 
in order to give more importance to the metabolites with fewer connections. The reason to do this is that 
metabolic networks of all organisms usually share their hubs (metabolites with many connections), so in 
order to establish differences and similitudes for different networks, one should focus on specific metabolites 
particular to only some organisms sharing common features. 

On the other hand, the factor |v^\Vb | ^-*v-ev A \v B ^ e s( v j) gives an average of the number of connections 
of the metabolites only present in A respect to the whole network. This is done in order to rescale the size 
of the network. 

Analogously, we can define /3 and 7 from the metabolites in the other two sets. 

\ Vj£V B \V A J Vi£V B \ A 

1= [wa^m £ de ^) E ggj 

For illustrating the process, let's consider three organisms, the Synechocystis sp. PCC 6803 ( that we refer 
to as syn), Synechococcus elongatus PCC7942 (referred to as syf), and the Escherichia coli K-12 MG1655 



(referred to as eco). In Table 5.4 one can see the number of metabolites and edges in the networks of these 
organisms, and in Table [53] we show the number of elements in each one of the three sets of the partition in 
which we split the set of vertex of the network obtained from each pair of these three organisms. 

Let's now focus on a few metabolites to see their contribution to the differentation parameters (i.e. to the 
parameters a, j3, and 7). For this, we chose Pyruvate (PYR), Glyoxylate (GXL), and 2-Dehydro-3-deoxy- 
6-phospho-D-gluconate (6PDG), which are respectively very, medium, and poorly connected metabolites 



present in these three organisms. In Table [5l3| we show the contribution of these metabolites to the parameters 
a, /?, and 7. 

Finally, the comparison between the networks A and B, namely (a,b, is defined as: 

, \v A \ a+ \v B \P 



20 

The parameters a and f3 are balanced since some organisms have much smaller metabolic networks than 
others. If this is not corrected, it results in a disproportionate size between subnetworks generated by Va\b 
and by Vb\a- In order to weaken this difference, the parameter factors 1^4 and jj^j are introduced. For 
two identical networks, a and /3 are zero, and so that £ = 0. For two networks which have not a single 
metabolite in common we have 7=0 and so £ = 00. 

3 Reconstruction of the phylogenetic tree 

Given a set of n organisms {Ax, A2, . . . , A n }, we will see how to construct their phylogenetic tree taking into 
account the degrees of similarity between every pair of metabolic models. 

Firstly, let TV = (V,E,w) be a non-directed connected complete weighted network, where V — 
{Ai, A2, ■ ■ ■ , A n } is the set of vertex that represent the metabolic models of the aforementioned organ- 
isms, E is the set of edges (Aj, Aj), 1 < i, j < n, i 7^ j, and w : E —> R is a function that assigns to every 
edge (Ai, Aj) the amount w^j — Ca^A - ■ Looking at the definition of Q, we observe that this network N must 
be symmetric. In particular, all the weights in our study are strictly positive. 

Secondly, we will compute a minimum spanning tree of N, that is a tree which has V as the set of vertex, 
and such that the sum of the weights associated to the edges of this tree is minimum. In these trees, every 
vertex Ai € V is connected with at least one of the other vertex of V \ {Ai} by an edge that has minimum 
weight among all the edges incident to Ai. The well-known Kruskal algorithm give us a procedure for finding 
these trees, see for instance [28]. We just have to follow the trace of Kruskal algorithm in order to recover 
the phylogenetic tree of the organisms represented by the models Ax,..., A n . 

In order to compute the phylogenetic tree of the models {Ax, A 2 , . . . , A n }, consider the minimum spanning 
tree of N, namely T = (V, E' , w\e<) where E' C E and w\e> denotes the restriction of the function w to the 
elements in E' . Let us take all the elements of E' in decreasing order of weights, that is E' = {e[, e' 2 , . . . , e' n _x} 
with w(e'x) > w(e' 2 ) ■ ■ ■ > w(e' 1 ). We are going to remove edges from T following this order. Everytime an 
edge is removed, the number of connected components of the resulting graph is increased in one respect to 
the previous one. We can represent this division of connected componets by a binary tree. The phylogenetic 
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tree is generated taking into account how we divide T. 

There are two different situations depending on the size of the (new) connected components (if any of 
them consists on a single vertex or not). Let us start with the edge with maximum weight in T, that we 
have denoted as e! x . Suppose that e' x is adjacent to two vertex Ai and Aj , with 1 < ig, jo < n, i ^ j. Then 
two situations can happen: 

(a) One of these vertex, for instance Aj , is a leaf (vertex of degree 1), 

(b) No one of these two vertex is a leaf (each vertex is still connected with other vertex). This happens 
only if the former connected component has 3 or more vertex. 

We point out that our phylogenetic tree will have two types of vertex: the leaves, that represent metabolic 
models, and the inner vertex, that represent that there are two branches and each one has more than one 
vertex. 

We start our phylogenetic tree with a vertex vq that will be its root. Then two vertex Vi 7 v 2 are hanged 
from v . Each one of these vertex represents one of the two connected components of the network T \ {e^}. 
Let us see what to do with v\ and v 2 depending in which case we are. 

• If we are in case (a) , one of these two vertex, for instance V\ , represents the vertex A io and v 2 the other 
one, i>2, represents the other connected component of T, which is generated subgraph of T generated 
by the vertex of V \ {A io }. 

• If we are in case (b), one of the vertex, for instance V\, represents the connected component of T\ {e^} 
that contains Ai , and the other vertex, v 2 , the connected component of T \ {e^} that contains Aj . 

This procedure is repeated again with v± and v 2 and removing e' 2 from T\ {e^}. When we remove e' 2 , 
then cither the connected component that represents v\ or v 2 is split into two smaller ones, and the vertex 
associated to this component plays again the role of vq. This proccess is repeated until we remove all the 
edges. 

Let us see how it works with two examples: 

1. In Table [5/7] we have the weights associated to a set of 10 organisms. The edge of maximum weight in 
the minimum spanning tree associated to the underlying network is the one that connects the mge with 
the lpl, with weight 0.123. We can see in Figure [572] that two vertex are hanging from the root of the 
tree. The one on the left represents the mge, the one on the right represents the subgraph associated 
to the rest of vertex, where the lpl can be found. 

G 



2. In the case of 33 organisms, when we remove from the minimum spanning tree the edge with maximum 
weight, we split this tree into 2 connected components: the one associated to the pair mge and mpm, 
and one the associated to the other vertex. 

Finally, the vertex in the phylogenetic tree can keep more information concerning the aforementioned 
minimum spanning tree. Suppose that the height of our phlyogenic tree is w(e^), that represents the 
maximum of the weights in the minimum spanning tree. i.e. the weight associated to e' Y - We place the root 
of our phylogenetic tree at height y — w(e[). Now, two vertex are hanged from the root. If one is associated 
to a single vertex for instance v\ in case (a), then we place this vertex at height y = 0. We remember that 
this vertex represents the metabolite Aj . If not, for instance v 2 in case (a) and either V\ or v 2 in case (b), 
each one of these vertex represents a connected component with more than one vertex in which the minimum 
spanning tree is split. In order to know at which height we should put these vertex we have to continue 
removing edges from the former tree. After removing e' 2 , one of these connected components, for instance 
the one represented by v 2l is split again into 2 smaller connected components. So we place the vertex v 2 at 
height w(e' 2 ). We repeat this process recursevely until the initial tree is just reduced to isolated vertex. 

4 Results 

We have reconstructed two phylogenetic trees, one with 10 bacteria and another one with both, prokaryotes 



and eukaryotes. In Table 5.7 we show the parameter £ for the pairwise comparison of the 10 prokaryotes in 
the first tree. The data for the comparison of the 33 organisms in the second tree is given in a file in the 
supplementary materials. 

The organisms in each comparison are: 

• 10 organisms tree — > Mycoplasma genitalium (mge), Lactobacillus plantarum WCFS1 (lpl), Syne- 
chocystis sp. PCC 6803 (syn), Synechococcus elongatus PCC7942 (syf), Synechococcus elongatus 
PCC6301 (syc), Clostridium beijerinckii (cbe), Burkhoderia cenocepacia J2315 (bcj), Escherichia coli 
K-12 MG1655 (eco), Thermotoga maritima (tma) and Yersinia pestis KIM10 (ypk). 

• 33 organisms tree — > Mycoplasma genitalium (mge), Mycoplasma pneumoniae 309 (mpm), Synechocys- 
tis sp. PCC 6803 (syn), Synechococcus elongatus PCC7942 (syf), Synechococcus elongatus PCC6301 
(syc), Clostridium beijerinckii (cbe), Salmonella bongori (sbg), Escherichia coli K-12 MG1655 (eco), 
Aquifex aeolicus (aae), Yersinia pestis KIM 10 (ypk), Cyanobacterium UCYN-A (cyu), Thermosyne- 
chococcus elongatus (tel), Microcystis aeruginosa (mar), Cyanothece sp. ATCC 51142 (cyt), Cyanoth- 
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ece sp. PCC 8801 (cyp), Gloeobacter violaceus (gvi), Anabaena sp. PCC7120 (ana), Anabaena azollae 
0708 (naz), Prochlorococcus marinus SS120 (pma), Trichodesmium erythraeum (ter), Acaryochloris ma- 
rina (amr), Halophilic archaeon (hah), Polymorphum gilvum (pgv), Micavibrio aeruginosavorus (mai), 
Agrobacterium radiobacter K84 (ara), Clostridiales genomosp. BVAB3 (clo), Gamma proteobacterium 
HdNl (gpb), Vibrio fischeri ESI 14 (vfi), Vibrio fischeri MJ11 (vfm), Haemophilus influenzae F3031 
(hif), Coprinopsis cinerea (cci), Sus scrofa (ssc) and Leishmania braziliensis (lbz). 



In Figures |5.2| and |5.3| we present the two phylogenetic trees that we have constructed. 

In the first tree the only organism displaced in relation to what is expected from standard methods of 
phylogenetic tree reconstruction is the tma. In both trees mge (and mpm in the second one) diverge from 
other organisms at the beginning of the tree. This happens because of their minimalistic genomes, with only 
a couple of hundred of metabolites in their metaboloms. As a result, when comparing with an organism 
without a reduced genome with around almost a thousand metabolites, several hundred metabolites will 
not have a correspondent one, increasing hugely the value of a in the calculation of the parameter £, and 
therefore distancing these organisms from the rest. In any case, we used for the second study organisms from 
very different origins in the evolutionary history and we found that the method is able to separate bacteria, 
archea and eukaryotes. Different strains of the same species also appear closely related and sharing branches 
with organisms from the same family and order. 

We have also studied the sensibility of the parameter £. For this we performed a monte-carlo analysis of 
£. The procedure for this analysis is explained as follows. Given two organisms, one of them remains the wild 
type while, with the other, one builds an ensemble with N t elements, where each element is the result of nx 
knock-outs (removal of tik randomly selected reactions from the metabolic model) in the organism. Then the 
calculation of Q is performed between the wild type organism and each organism in the knock-out ensemble. 
From this process one obtains an ensemble of Nt values of £ for the comparison (one from each version of the 
organism in the knock-out ensemble), from which one calculates its average and standard deviation. This 
standard deviation is treated as an indicator of the sensibility of the parameter (as a function the number 
of knock-outs). 

We performed this sensibility analysis for four organisms (syn, syf, eco and mge) with ensembles of sizes 
N t — 500 for uk =5, 10, 50 and 100. The results are in Tables |5.8||5.11 These four organisms have been 



chosen in order to observe the sensibility in the comparison between very similar organisms (syn and syf), 
more distant ones (syn and eco) and very different ones (syn and mge). 
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5 Conclusions 

In this work we have developed a methodology for comparing organisms based on their metabolic networks. 
This methodology has been successfully applied for the reconstruction of phylogenic trees for several organ- 
isms from a broad range of families and kingdoms. Resulting trees stand well their comparison with the 
so-called "tree of life" . The great majority of the branches in the tree fit well their expected positions and 
their distance is in good correlation with evolutionary distances. The discrepancies found can be explained 
by particularities in these very few organisms not fitting the tree, such as tremendous genome reductions 
that caused reduced metabolisms. 

Our methodology is innovative for it is not directly based on the structure and evolution of proteins or 
DNA, but on the metabolism and on the organisms' components and metabolic capabilities, allowing one 
to compare organisms very distant from the evolutionary point of view or organisms for which ortholog's 
comparison is difficult. In order to accomplish this, we make use of the correlation between evolutionary 
distances and metabolic network likelihood and propose our methodology as a starting point to study it. 

Metabolism information is retrieved as a subset of the whole genome information. We hereby show that 
metabolic network connectivity can be used to build phylogenetic trees that are in accordance with gene- 
directed trees. It can be argued whether the selected construction parameter (£) is the optimal one for this 
purpose (or even if there is an optimal one), but it stands clear that this is an innovative application for 
metabolic models, their curation and cross-species evolutionary studies. 
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Figures 

5.1 Figure 1 - conjs.pdf 
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(Color on-line) Representation of the sets of metabolites between two organisms. 
5.2 Figure 2 - treelO.png 
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mge syn syf syc cbe Ipl tma bcj eco ypk 



(Color on-line) Phylogenetic tree with 10 organisms. 
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5.3 Figure 3 - tree30.png 



n 



(Color on-line) Phylogenetic tree with 33 organisms. 



Tables 

5.4 Table 1 - Vertex and edges in the networks of syn, syf and eco. 



Organism 


# Vertex 


# Edges 


syn 


1001 


2891 


syf 


979 


2810 


eco 


1227 


3801 



5.5 Table 2 - Metabolites in the three sets of the partition when comparing three organisms. 





syf 


eco 




\v A nv B \ 


= 911 


\v A n v B 


= 778 


syn 


\Va\V b \ 


= 90 


\Va\Vb\ 


= 223 




\Vb\V a \ 


= 68 


\Vb\Va\ 


= 449 


syf 




\v a n v B 

\Va\V b \ 
\Vb\Va\ 


= 775 
= 204 
= 452 
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5.6 Table 3 - Contribution of different metabolites to the differentiation parameter (() between two 
networks. The column shows the weight of the metabolite in the calculation of PAnB.i which is 
the inverse of the degree of the metabolite divided by the sum of the inverses of the degrees of 
all metabolites contributing to the parameter. 



Metabolite 


Organisms in 
comparison 


PAnB.i 


6i 


Contribution (%) 




syn and syf 


0.98 


0.127 


0.0064 


PYR 


syn and eco 


0.73 


0.117 


0.0044 




syf and eco 


0.75 


0.113 


0.0044 




syn and syf 


0.86 


0.454 


0.020 


GXL 


syn and eco 


0.87 


0.550 


0.024 




syf and eco 


0.80 


0.439 


0.018 




syn and syf 


1.00 


3.176 


0.16 


6PDG 


syn and eco 


0.80 


1.762 


0.072 




syf and eco 


0.80 


1.757 


0.072 



5.7 Table 4 - Comparison matrix for ten organisms. 



org 


mge 


lpl 


syn 


syf 


syc 


cbe 


bcj 


eco 


tma 


ypk 


mge 


0.0 


0.123 


0.1951 


0.1753 


0.1906 


0.1438 


0.1374 


0.1384 


0.1306 


0.1384 


lpl 


0.123 


0.0 


0.1554 


0.1543 


0.1645 


0.0719 


0.1217 


0.1227 


0.0716 


0.1154 


syn 


0.1951 


0.1554 


0.0 


0.0195 


0.0188 


0.1286 


0.1135 


0.1137 


0.1687 


0.1248 


syf 


0.1753 


0.1543 


0.0195 


0.0 


0.0054 


0.1262 


0.117 


0.1076 


0.1677 


0.1174 


syc 


0.1906 


0.1645 


0.0188 


0.0054 


0.0 


0.1327 


0.1149 


0.1043 


0.1657 


0.1133 


cbe 


0.1438 


0.0719 


0.1286 


0.1262 


0.1327 


0.0 


0.1121 


0.0933 


0.0737 


0.1071 


bcj 


0.1374 


0.1217 


0.1135 


0.117 


0.1149 


0.1121 


0.0 


0.0678 


0.1324 


0.068 


eco 


0.1384 


0.1227 


0.1137 


0.1076 


0.1043 


0.0933 


0.0678 


0.0 


0.1101 


0.0295 


tma 


0.1306 


0.0716 


0.1687 


0.1677 


0.1657 


0.0737 


0.1324 


0.1101 


0.0 


0.1075 


ypk 


0.1384 


0.1154 


0.1248 


0.1174 


0.1133 


0.1071 


0.068 


0.0295 


0.1075 


0.0 



5.8 Table 5 - Sensibility calculation for N t = 500 and % = 5. Each element in the table is the average 
of the parameter ( in an ensemble plus (minus) its standard deviation ((±a^). 



org / org 


syn 


syf 


eco 


mge 


syn 


0.0002 ± 0.0003 


0.0184 ± 0.0005 


0.0893 ± 0.0005 


0.1600 ± 0.0014 


syf 


0.0184 ± 0.0004 


0.0002 ± 0.0003 


0.0857 ± 0.0006 


0.1527 ± 0.0014 


eco 


0.0892 ± 0.0005 


0.0856 ± 0.0005 


0.0001 ± 0.0002 


0.1278 ± 0.0009 


mge 


0.1597 ± 0.0025 


0.1527 ± 0.0026 


0.1283 ± 0.0015 


0.0014 ± 0.0016 



5.9 Table 6 - Sensibility calculation for N t — 500 and % = 10. Each element in the table is the 
average of the parameter ( in an ensemble plus (minus) its standard deviation ((±a^). 



org / org 


syn 


syf 


eco 


mge 


syn 


0.0005 ± 0.0005 


0.0186 ± 0.0006 


0.0896 ± 0.0007 


0.1604 ± 0.0018 


syf 


0.0187 ± 0.0006 


0.0005 ± 0.0005 


0.0860 ± 0.0007 


0.1532 ± 0.0019 


eco 


0.0893 ± 0.0008 


0.0857 ± 0.0007 


0.0003 ± 0.0003 


0.1281 ± 0.0011 


mge 


0.1602 ± 0.0035 


0.1531 ± 0.0032 


0.1288 ± 0.0023 


0.0028 ± 0.0023 
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5.10 Table 7 - Sensibility calculation for N t = 500 and % = 50. Each element in the table is the 
average of the parameter ( in an ensemble plus (minus) its standard deviation ((±a^). 



org / org 


syn 


syf 


eco 


mge 


syn 


0.0028 ± 0.0011 


0.0209 ± 0.0014 


0.0915 ± 0.0017 


0.1652 ± 0.0045 


syf 


0.0207 ± 0.0013 


0.0029 ± 0.0011 


0.0879 ± 0.0016 


0.1575 ± 0.0044 


eco 


0.0903 ± 0.0017 


0.0868 ± 0.0016 


0.0016 ± 0.0007 


0.1301 ± 0.0029 


mge 


0.1638 ± 0.0080 


0.1577 ± 0.0077 


0.1343 ± 0.0055 


0.0170 ± 0.0053 



5.11 Table 8 - Sensibility calculation for N t = 500 and hk = 100. Each element in the table is the 
average of the parameter ( in an ensemble plus (minus) its standard deviation ((±a^). 



org / org 


syn 


syf 


eco 


mg 


e 


syn 


0.0058 ± 0.0016 


0.0239 ± 0.0020 


0.0942 ± 0.0024 


0.1715 


± 


0.0062 


syf 


0.0238 ± 0.0018 


0.0061 ± 0.0017 


0.0907 ± 0.0023 


0.1630 


± 


0.0066 


eco 


0.0919 ± 0.0024 


0.0883 ± 0.0022 


0.0033 ± 0.0011 


0.1329 


± 


0.0040 


mge 


0.1694 ± 0.0120 


0.1648 ± 0.0131 


0.1433 ± 0.0092 


0.0460 


± 


0.0076 



Additional Files 

Additional file 1 - comps_33orgs.dat 

Table with the parameter £ resulting from the comparison of 33 organisms. 
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